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Abstract 

We investigate the possibility of superfluid/ferromagnet/superfluid (SFS)-junction in a superfluid 
Fermi gas. To examine this possibility in a simple manner, we consider an attractive Hubbard 
model at T = within the mean-field theory. When a potential barrier is embedded in a superfluid 
Fermi gas with population imbalance (iV-j- > JVj,, where N a is the number of atoms with pseudospin 
a =t, i), this barrier is shown to be magnetized in the sense that excess t-spin atoms are localized 
around it. The resulting superfluid Fermi gas is spatially divided into two by this ferromagnet, so 
that one obtains a junction similar to the superconductor/ferromagnet/superconductor-junction 
discussed in superconductivity. Indeed, we show that the so-called 7r-phase, which is a typical 
phenomenon in the SFS-junction, is realized, where the superfluid order parameter changes its 
sign across the junction. Our results would be useful for the study of magnetic effects on fermion 
superfluidity using an ultracold Fermi gas. 

PACS numbers: 03.75.Ss, 03.75. +b, 03.70. +k 
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I. INTRODUCTION 



The advantage of recently realized superfluid 40 K and 6 Li Fermi gases 1-4] is that one 
can experimentally tune various physical parameters. This unique property enables us to 
study various physical properties of Fermi superfluids in a wide parameter region^. Since 
the background physics of superfluid Fermi gas is essentially the same as superconductiv- 
ity, the study of former atomic system is also useful for understanding the latter electron 
system. Indeed, using a tunable interaction associated with a Feshbach resonance[l-4, 6], 
one can study superfluid properties from the weak-coupling regime to the strong-coupling 
limit in a unified manner [7- 12]. In the intermediate coupling re gim e (crossover region), 
strong pairing fluctuations lead to the pseudo gap p henomenon 13Hl7j , which is related to 



the under-doped regime of high-T c cuprates 



18H24J. In addition, a cold Fermi gas loaded 

which is one of the 



25 



onto an optical lattice is well described by the Hubbard model 
most fundamental models in strongly-correlated electron system. Thus, using this lattice 
Fermi gas, we may study various phenomena originating from strong correlation, such as 



Mott insulator 



291 ]. and unconventional pairing state 30|. We also point out that a superfluid 



Fermi gas with population imbalance is similar to superconductivity under an external mag- 
netic field, so that the possibility of Fulde-Ferrel-Larkin-Ovchinnikov(FFLO) state 3lM33] 



34 ^ 



has been recently discussed in polarized Fermi gases 

The above examples clearly show the importance of cold Fermi g clS cLS db quantum sim- 
ulator for metallic superconductivity. However, at the same time, we also know that the 
current stage of superfluid Fermi gas cannot cover all the important topics discussed in 
superconductivity. In metallic superconductivity, magnetism is known to strongly affect 
superconducting properties 39l. |40|. In particular, the conventional s-wave pairing is easily 



destroyed by ferromagnetism, because the latter tends to align spins in the same direction, 

39(. On the other hand, in superfluid 



which is unfavourable for the singlet Cooper pairing 
Fermi gases, 'spins' in Cooper pairs are actually pseudospins describing two atomic hyper- 
fine states. Thus, at present, there is no object corresponding to the magnetic impurity in 
superconductivity. Thus, if we can introduce a magnetic object to a cold Fermi gas, together 
with the tunable interaction and optical lattice, this system would become a more useful 
quantum simulator to study metallic superconductivity. 

In this paper, we theoretically discuss an idea which enables us to study magnetic effects 
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on superfluid Fermi gases. Among various magnetic effects discussed in metallic supercon- 
ductivity, in this paper, we focus on a superffuid/ferromagnet /superfluid (SFS)-junction. 
In superconductivity, the SFS-junction in FigfTJ^a) is well known to modulate the phase of 



superconducting order parameter. Namely, it induces the 7r-phase[4lH45| shown in Figfjjb), 
where the order parameter changes its sign across the SFS-junction. Such a phase mod- 
ulation does not occur in the superconductor/insulator/superconductor (SlS)-junction in 
Figfjjc), where what we call the 0-phase shown in Figfjjd) is realized. 

The origin of 7r-phase in the superconducting SFS-junction is deeply related to the FFLO 
state. When the order parameter A penetrates into the ferromagnetic junction, a molecular 
field in the ferromagnet induces an FFLO-like spatial oscillation of A. Inside the junction, 
when the sign of oscillating order parameter at the right edge of the junction is opposite to 
the sign of the order parameter at the left edge, the 7r-phase is realized. On the other hand, 
when the order parameter has the same sign at both the edges, the 0-phase is obtained. 
Indeed, a transition between 0-phase and 7r-phase (0-7T transition) has been observed in 
Nb/ Cun.47Nio.53/Nb-j unction, when one varies the thickness of ferromagnetic Cuo.47Nio.53- 
layer 45 ] . 

We explain the outline of our idea. In a superfluid Fermi gas with population imbalance, 
the phase separation has been observed 46|, |47|, as schematically shown in Figj2]^a). Then, 
noting that excess atoms tend to occupy the spatial region where potential energy is large 



541 ]. we expect that, in a double well potential, excess atoms are also localized around the 



central potential barrier, as shown in Figj2]^b). Regarding this central region as a ferromagnet 
of pseudospin-'f, Figf2](b) looks like an SFS-junction. Although this ferromagnet is actually 
a gas of atoms in the hyperfine state described by pseudospin-7, we will find that it indeed 
induces the rr-phase as in the case of real ferromagnet. 

This paper is organized as follows: In Sec. II, we explain our formulation. Treating an 
attractive Hubbard model within the mean field theory at T = 0, we numerically determine 
the superfluid order parameter, as well as the particle density, around a potential barrier 
in the presence of population imbalance. We compare the energy of 7r-phase with that of 
0-phase to examine the stability of 7r-phase. In Sec. Ill, we examine the possibility of SFS- 
junction and 7r-phase in a simple one-dimensional model in the absence of a trap. Effects of 
a harmonic trap, as well as dimensionality, will be discussed in Sec. IV. 
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FIG. 1: (a) Schematic SFS-junction. 'S' and 'F' mean the superfluid (superconducting) region 
and ferromagnet, respectively, (b) Spatial variation of superfluid order parameter A(x) in the 
7r-phase. (c) SIS-junction, where T means the insulating barrier, (d) Spatial variation of A(x) in 
the SIS-junction (0-phase). 



II. FORMULATION 

We consider a two-component Fermi gas with population imbalance. To examine the 
possibility of SFS-junction and 7r-phase in a simple manner, we employ the Hubbard model 
described by the Hamiltonian, 

H = -tJ2 + h - c - - U Yl "*.t"U + ^2i V i- hi,*- (!) 

(i,j),cr i i,<r 

Here, cj )(T is the creation operator of a Fermi atom at the i-th site, where the pseudospin 
a =t, | describes two atomic hyperfine states, t is a nearest-neighbour hopping, and the 
summation is taken over the nearest-neighbour pairs. —U{< 0) is an on-site pairing 
interaction, and hi )(T = c ia Ci t(r is the number operator at the i-th site. Since we consider 
a polarized Fermi gas, the chemical potential fi a depends on pseudospin a. Vi is a (non- 
magnetic) potential, consisting of a potential barrier V t b and a harmonic trap V*. Their 
detailed expressions are given later. 

Although a real superfluid Fermi gas is a three-dimensional continuum system, we con- 
sider a one-dimensional chain and two-dimensional square lattice in this paper. In this 
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FIG. 2: (a) Schematic density profile of a trapped superfluid Fermi gas with population imbalance 
(iV-f > iVjJ. Excess t-spin atoms are localized at the edges of superfluid region (ti)- (b) Density 
profile expected in a double well potential. Excess atoms can be also localized around the central 
potential barrier. When we regard the central spin-polarized region as a ferromagnet, the system 
can be viewed as an SFS-junction. 



regard, we emphasize that the assumed low dimensionality and lattice structure are not es- 
sential for the realization of SFS-junction and 7r-phase discussed in this paper. However, to 
suppress band effects by the lattice, we consider the case of low particle density (n iiCT ) < 0.3. 

We treat Eq.([T]) within the mean field theory at T = 0. Under this approximation, we 
have 



H 



MF 



^ [Vj - l*>a - U (hi -a)] n i>a + 



u 



+ [/(n i;t )(n a ) 



(2) 



Here, A, = U (ci^Ci^) is the superfluid order parameter (which is taken to be real in the 
following calculations). Equation (TSJ) can be rewritten as 



H MF = + 



A 2 

Vi-m-U (n iit ) + + U (n iit ) (h lA ) 



(3) 



Here, = [c\ f c|,t 



• • • c M + 624 • • • cm,J, where M is the total number of lattice sites. 
The matrix elements of H are chosen so that Eq.([3]) can reproduce Eq.(j2]). As usual, Eq.((3]) 



5 



can be conveniently diagonalized by the Bogoliubov transformation, 



/ 6l , \ 

C2,f 



CM,t 
C f 

C ] 
2,4. 



a 2 ,t 



t 



a 



24 



(4) 



V fi M4 / V & M4 / 

where W is a 2M x 2M-orthogonal matrix. The diagonalized Hamiltonian has the form 



M 



3, a 



A 2 



(5) 



where _Ej j0 - is the Bogoliubov excitation spectrum. The superfluid order parameter Aj 
U(ci^Ci^) and particle density (n i)0 .) in Eq.(J5]) are given by, respectively, 



M 

,M+j WM+i,M+j8(Ej,i) 



(6) 



M 



A/ 



(7) 



In the presence of population imbalance, since Ej >a may be negative, the step functions 
appear in Eqs.flH]) and ([7]). In the unpolarized case, one can safely take Ej^ = Ej^ > [39 ] . 

We numerically diagonalize H, and self-consistently determine A i; (n it(r ), as well as fi a , 
for a given parameter set (U, N^, N±, Vi), where N a = Yli{™i,o) * s the total number of cx-spin 
atoms. Once the self-consistent 7r-phase solution is obtained, we energetically compare this 
state with the 0-phase solution. In this procedure, since these two self-consistent solutions 
usually have different values of chemical potentials fi a , it is convenient to use the free energy 
Eq rather than the thermodynamic potential Q (= (Hmf) at T = 0). The former is imme- 
diately obtained from the latter by the Legendre transformation, as Eq = Q + ^2 a fi a N a . 
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At T = 0, we have 



Vi-ni-U (h i7t ) + JJ- + U (n iit > (n a > - E hi 



(8) 



We will use the energy difference AEq = ^ _phafaC — ^o_-phase ] 3e ^ ween ^he -^.phase and 
0-phase to examine the stability of 7r-phase. 

Here, we present the detailed expression for the potential Vt = V* + Vf. In Sec. Ill, we 
ignore effects of a trap V? for simplicity, and only take into account the potential barrier in 
the one-dimensional system. However, considering the model with only one potential barrier 
under the periodic boundary condition, we cannot obtain the 7r-phase solution, because it 



does not satisfy the periodic boundary condition imposed at the edges of the system. To 
avoid this problem, in this paper, we consider the system with two potential barriers shown 
in FigJH^a). Namely, we set 

( V b e -(i x -L*/Af/e \<i x< LJ2, 
\ Vftr^-aW*)"/* L x /2 <i x < L x , 

where L x is the number of lattice sites in the x-direction, and i x is the site number. (The 
lattice constant is taken to be unity.) Vq and I describe the height and width of the potential 
barrier, respectively. We briefly note that the assumed Gaussian spatial dependence in 
Eq.Q is not crucial for the formation of SFS-junction and 7r-phase. Using Eq.Q, we 
can conveniently obtain both the 0-phase and 7r-phase solutions under the same periodic 
boundary condition. We also use this prescription in considering the two-dimensional model 
in FigE(b). 

In Sec. IV, we include the harmonic trap V*. In the one- dimensional chain, we set V, as 

Vi = VS (i*-^) 2 + V b e~^- L *W e , (10) 



2 , 

where Vq describes the strength of harmonic trap. The potential barrier is centered at 
i x = L x /2. In the case of two-dimensional cigar trap, we take 



V 



-Li x \ „ / . Li. 



+ v b e -(i*-L*m>/e^ 

where the cigar trap is described by (Vq x , Vq v ), and L y is the number of lattice sites in the 
y-direction. In Eq. flTT]) . the potential barrier V^ is centered at i x = L x /2 and parallel to the 
y- axis. 



III. FORMATION OF SFS-JUNCTION AND vr-PHASE 

Figure H] shows the self-consistent 7r-phase solution (a) and 0-phase solution (b) in the 
one- dimensional model in Fig|3^a). As expected, the polarization Sj = (n^) — (^4) is 
finite around the barrier, which may be viewed as a ferromagnetic junction. In FigJH (where 

— Ni = 2), since there are two potential barriers in the model, the number of excess 
atoms localized around each barrier (= Nl) equals one. Strictly speaking, there is no clear 
boundary between the superfluid region and spin-polarized region because of the continuous 
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FIG. 4: Self-consistent 7r-phase solution (a) and 0-phase solution (b) in the one-dimensional model 
in FigG^a). From the comparison of their energies, we find that the 7r-phase is more stable than 
the 0-phase. We also show the case of unpolarized Fermi superfluid in (c). The particle density 
(nj j(J ) and superfluid order parameter Aj are shown in the upper and middle panels, respectively. 
For convenience, the potential barrier is also shown in the lower panels. The total number of 
atoms equals 160, and we take Nf = 8L_ N± = 79 for (a) and (b), and Nf = N i = 80 for (c). The 
other parameters are taken as U/t = 4(55|], L x = 300, Vq ft = 0.1 and I = 5. The same system size 
and parameters for the barrier are also used in FigsJMS Since results are always symmetric with 
respect to i x = L x /2, we only show the region 1 < i x < L x /2 here. 
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variations of Aj and s«. In this paper, we simply call the region where Si is large to some 
extent the ferromagnetic junction. For example, in the case of FigJU from the spatial 
variation of Sj shown in FigJSl we roughly regard the region 72 < i x < 78 as the ferromagnetic 
junction. 

In the 7r-phase solution (FigJU^a)), the order parameter Aj changes its sign across the 
ferromagnetic junction, which is in contrast to the 0-phase case (FigJU^b)). From the com- 
parison of their energies, the 7r-phase is found to be more stable than the 0-phase. (Note 
that AEq < at Nl = 1 in Figj6j) This confirms that the spin-polarized region really works 
as a ferromagnetic junction. 
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FIG. 5: Spatial variation of polarization Sj = — (n^j.) around the barrier, (a) shows effects of 
the number of localized excess atoms Nl- (b) shows effects of the pairing interaction U. 'V ("0") 
in the legend means that the 7r-phase (0-phase) is the the most stable state in each case. 

We note that the 7r-phase is not obtained in the absence of population imbalance. As an 
example, FigJU^c) shows the case of N-f = N± = 80, where the order parameter (panel (c2)) 
is found to be not affected by the barrier, except for a slight decrease around it due to the 
slight decrease of particle density shown in FigH^cl) js6|. 

The reason why excess ^-spin atoms are localized around the barrier is understood as 
follows. When the phase separation occurs in a polarized Fermi superfluid, the superfluid 
order parameter is destroyed in the spatial region where excess atoms are localized. Then, 
to make the loss of condensation energy by this pair-breaking effect as small as possible, the 
localization should occur around the potential barrier, where the order parameter is already 
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FIG. 6: Energy difference AEq = phase _ phase between the 7r-phase and 0-phase, as a 
function of the number of localized excess atoms per one potential barrier. Since there are two 
barriers in FigEfa), N L equals [iY t - iVj/2. We take U/t = 4 and N = 160. 



(slightly) suppressed in the unpolarized shown in FigH^c2). 

As explained in the introduction, the 7r-phase originates from the FFLO-type spatial 
oscillation of order parameter in the ferromagnetic junction. Because of this mechanism, 
the 0-7r transition has been observed in superconducting Nb / Cun 47NI0 ^ /Nb-j unction |45] . as 
one varies the thickness of ferromagnetic Cuo.47Ni .53-layer. In the present superfluid Fermi 
gas, since the thickness of ferromagnetic junction depends on N L (See FigfS^a).), we can also 
expect the same phenomenon. Indeed, Figj6] shows that the 0-phase becomes more stable 
than the 7r-phase when two excess 'f-spin atoms are localized around the potential barrier 
(Nl = 2). We show the self-consistent solutions in this case in FigJTl In this case, since the 
0-phase is the most stable state, to make the 7i~-phase solution, one needs to forcibly bend 
the order parameter in the ferromagnetic junction, which leads to the fine structure around 
the barrier seen in FigJ7]^a2). 

Figure [6] shows that the 7r-phase again becomes more stable than the 0-phase when 
Nl = 3. In this case, we see in FigJHJ^b) an FFLO-like oscillation with the wavelength A ~ 9 
in the ferromagnetic junction (70 < i x < 80). In the ordinary FFLO state, the wavelength 
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FIG. 7: (a) 7r-phase and (b) 0-phase solutions when Nf = 82, N± = 78, and U/t = 4. In this case, 
two excess atoms are localized around the potential barrier, and the 0-phase is more stable than 
the 7r-phase. 



of the order parameter oscillation is given by 33] , 

v F 



A ~ lu- 



ll 



(12) 



where vp and h are the Fermi velocity and an external magnetic field, respectively (where 
h involves the Bohr magneton). In the present model Hamiltonian in Eq.(J2]), h in the 
ferromagnetic junction is evaluated as, 



(13) 



where s7 is the polarization averaged over the junction. Substituting numerical values into 
Eqs.( TT2|) and (|T3|) we obtain A ~ 9.1, which is consistent with the value (A ~ 9) evaluated 
from the inset in FigJH^b). 

When the ferromagnetic junction becomes very thick, the Josephson coupling between 
the two superfluid regions becomes very weak. Since the energy difference AEq between the 
7r-phase and 0-phase purely comes from the Josephson coupling, |A.E(j| in Fig 16] becomes 
small with increasing Nl. 
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FIG. 8: Stable vr-phase solution when N L = 3. We take iV t = 83, N i = 77, and U/t = 4. The inset 
shows the magnified spatial variation of Aj in the ferromagnetic junction (70 <, i x <, 80). 
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FIG. 9: Energy difference AEq between 7r-phase and 0-phase as a function of U. We take = 81 
and iVj_ = 79. The inset shows AEq magnified around the 0-tt transition. 



While the thickness of the ferromagnetic junction depends on Nl, Figj5]^b) shows that it 
is not so sensitive to the pairing interaction U. However, even in this case, Figj9] shows that 
the 0-7T transition occurs when U becomes large. This is because, when the BCS excitation 
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FIG. 10: Energy difference AEq as a function of Nl/L v , in the two-dimensional model. The x-axis 
means the number of localized excess atoms around the barrier per unit length in the y-direction. 
We set U/t = 7, V^/t = 0.2, £ = 1, and L x = 50, which are also used in Figjlll For the number of 
atoms, we take N a = (N + 2aN]f)/2 and N = 6L y . The value of L y (= 17 ~ 48) is chosen so that 
L v /Nl can be an integer. 

gap becomes large due to strong pairing interaction, the effective magnetic field given by 
Eq.( )13p also becomes large to realize a finite population imbalance [581] . Thus, the FFLO 
wavelength A in Eq.( Tl"2]) becomes short, leading to the transition from the 7r-phase to the 0- 
phase. We briefly note that, while the 0-ir transition occurs by tuning Nl and U, it does not 
occur when the potential shape is varied by tuning Vq and £, at least within our numerical 
calculations (although we do not show the results here). This means that the dominant 
role of potential barrier is a pinning center. Once excess atoms are pinned at the potential 
barrier, the magnetic effect is dominated by the number of localized atoms Nl, as well as 
the pairing interaction U. 

IV. EFFECTS OF DIMENSIONALITY AND TRAP POTENTIAL 

In the previous section, we considered the one-dimensional model in the absence of a 
trap, to simply examine the possibility of SFS-j unction and 7r-phase in a superfluid Fermi 
gas. In this section, we extend the previous study to include higher dimension and harmonic 
trap. 
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FIG. 11: (a) Stable O-phase solution when Ni/L y = 0.1 (L y = 30). (b) Stable 7r-phase solution 
when Ni/Ly = 0.25. (L y = 20). The upper and lower panels show the spatial variations of st and 
Aj, respectively. We only show the spatial region 1 < i x < 25 and 1 < i y < 20 in this figure. 



A. SFS-junction and 7r-phase in two dimension 

In the one-dimensional model discussed in Sec. Ill, although the SFS-junction is obtained, 
the stability of 7r-phase is sensitive to the number Nl of excess atoms localized around the 
potential barrier. This means that we have to adjust Nl within the accuracy of 0(1) 
in order to realize the 7r-phase, which is actually difficult experimentally. However, in a 
higher dimensional system, since excess atoms can align in the y- and ^-direction without 
remarkably increasing the thickness of ferromagnetic junction, the 7r-phase is expected to 
be more accessible. 

To see this, we examine the two-dimensional model shown in FigJH^b). As expected, we 
obtain the finite region of stable 7r-phase in FigfTn] (N L /L y > 0.12). Comparing FigsJTlTal) 
with (bl), we find that, while the increase of localized atoms results in the enhancement 
of polarization s« at the center of the barrier, the thickness of the ferromagnetic region is 
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FIG. 12: One-dimensional polarized superfluid Fermi gas trapped in a harmonic potential. The 
height of the potential barrier V b is taken to be (a) 0.4i, (b) 0.5i, and (c) 0.6i. We take = 47, 
iV|_ = 43, U/t = 2, Vq = 0.00005t, and i = 12. The density profile (nj j(J ), polarization S{, and A< 
are shown in the upper, middle, and lower panels, respectively. The potential barrier is centered 
at i x = L x /2 = 150. 



almost unchanged. However, since the junction would eventually become thick when one 
increases the number of excess atoms very much, the upper bound of the 7r-phase region in 
FigJTUI would exist at a large value of N L /L y (> 0.25). 



B. Effects of harmonic trap 

In a trapped system, excess atoms are not always localized around the potential barrier. 
Indeed, as shown in FigJT2l where the trap potential in Eq.f llOp is used, excess atoms are only 
localized at the edges of the gas cloud when the central barrier is low (panels (al) and (a2)). 
In this case, the potential barrier is non-magnetic, and the 0-phase is realized, as shown in 
panel (a3). As one increases the barrier height, the barrier is magnetized in the sense that 
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FIG. 13: (Color online) Two-dimensional polarized superfluid Fermi gas in a two-dimensional cigar 
trap. The central barrier height Vq is given by (a) 0.6t, (b) l.Ot, and (c) 1.2t. For other parameters, 
we take N t = 32, N± = 28, U/t = 4, V^Jt = 0.001, V^/t = 0.1, and I = 6. The upper panels 
show the total density profile, and the middle panels show the polarization. The order parameter is 
shown in the lower panels. In panels (b), the number of localized excess atoms around the central 
barrier equals one, and 7r-phase is obtained as the most stable state. In panels (c), the number of 
excess atoms around the barrier equals two, where the 0-phase is realized. For clarity, we explicitly 
write the sign of the order parameter ("+" and "-") in the lower panels. 



an excess atoms is localized around it, leading to the 7r-phase (panels (b)). The number 
of localized atoms around the central barrier further increases, as one further increases the 
potential height. In panels (c), the 0-phase is again realized due to the thick ferromagnetic 
junction by two excess atoms. 

The competition between the localization of excess atoms at the edges of the trap and 
magnetization of the central potential barrier can be also seen in the two-dimensional cigar 
trap given by Eq. fllip . We show the result for this case in FigJTHl Since the present system 
size is rather small due to computational problems, the 7r-phase in panels (b) soon changes 
into the 0-phase shown in panels (c), when the number of excess atoms localized around the 
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FIG. 14: Phase diagram of trapped polarized Fermi superfluid in terms of potential height Vq 
and width I. (a) One-dimensional harmonic trap, (b) Two-dimensional cigar trap. In this figure, 
'normal junction' is the region where excess atoms do not form a ferromagnetic junction in the 
trap center. Parameters in each case are the same as those in Figsll2l and 1131 respectively. 

central barrier increases only by one. However, as discussed in Sec. IV. A, this sensitivity 
would be suppressed when we consider a lager system. 

Figure HH shows the phase diagram in terms of the potential height and width. In both 
the one- and two-dimensional cases, we obtain finite regions of 7r-phase. Thus, the 7r-phase 
is expected to be experimentally accessible even for a trapped superfluid Fermi gas. 

V. SUMMARY 

To summarize, we have investigated the possibility of SFS-junction in a superfluid Fermi 
gas. Using the phase separation of a polarized Fermi gas, we showed that a non-magnetic 
potential barrier embedded in the system can be magnetized in the sense that some of excess 
atoms are localized around it. The polarized barrier works as a ferromagnetic junction, 
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leading to the 7r-phase. We also discussed how the presence of harmonic trap, as well as the 
dimensionality of system, affect the SFS-junction and 7r-phase. 

For the observation of 7r-phase, since the phase of the order parameter differs by tt across 
the junction, an interference experiment might be useful. Since it is known that low-energy 



bound states may appear around the nodes of order parameter [59j, |60|, the observation of 
these ingap states in single-particle excitation spectra may be another idea to detect the 
7r-phase. 

Besides the importance of the 7r-phase itself, it is also useful for various applications. For 
example, since it may be viewed as a part of FFLO state, one might use it to realize the 



FFLO state in cold Fermi gases. In addition, when the SFS-junction and 7r-phase are rea 



61 



ized 



62 ] 



in a superffuid Fermi gas confined in a toroidal trap, the so-called 7r-junction SQUID 
may be realized in a superfluid Fermi gas. Although we need further careful analyses for 
them, we expect that our results would be useful for the study of magnetic effects on Fermi 
superfluids by using cold Fermi gases. 
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